Low temperature ferromagnetism in (Ga, Mn)N 
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The magnetic properties of dilute magnetic semiconductors (DMSs) are calculated from first- 
principles by mapping the ab initio results on a classical Heisenberg model. By using the Korringa- 
Kohn-Rostoker coherent potential approximation (KKR-CPA) method within the local density ap- 
proximation, the electronic structure of (Ga, Mn)N and (Ga, Mn)As is calculated . Effective ex- 
change coupling constants Jij's are determined by embedding two Mn impurities at sites i and j in 
the CPA medium and using the Jij formula of Liechtenstein et al. It is found that the range of the 
exchange interaction in (Ga, Mn)N, being dominated by the double exchange mechanism, is very 
short ranged due to the exponential decay of the impurity wave function in the gap. On the other 
hand, in (Ga, Mn)As where p-d exchange mechanism dominates, the interaction range is weaker 
but long ranged because the extended valence hole states mediate the ferromagnetic interaction. 
Curie temperatures (Ic's) of DMSs are calculated by using the mean field approximation(MFA), 
the random phase approximation(RPA) and the, in principle exact, Monte Carlo method. It is 
found that the Tc values of (Ga, Mn)N are very low since, due to the short ranged interaction, 
percolation of the ferromagnetic coupling is difficult to achieve for small concentrations. The MFA 
strongly overestimates Tc- Even in (Ga, Mn)As, where the exchange interaction is longer ranged, 
the percolation effect is still important and the MFA overestimates Tc by about 50-100%. 



Dilute magnetic semiconductors (DMSs), such as (In, 
Mn)As and (Ga, Mn)As discovered by Munekata et al. 
and Ohno et al., have been well investigated as hopeful 
materials for spintronicsi. Curie temperatures (Ic's) of 
these DMSs are well establishediiSi^ and some prototypes 
of spintronics devices have been produced based on these 
DMSs. The magnetism in these DMSs are theoretically 
investigated and it is known that the ferromagnetism in 
these systems as well as (Ga, Mn)Sb can be well described 
by Zener's p-d exchange interaction, due to the fact that 
the majority d-states lies energetically in the lower part 
of the valence band^. Dietl et al^ and MacDonald et 
al.^ successfully explained many physical properties of 
(Ga, Mn)As based on the p-d exchange model, and first- 
principles calculations by Sato et al. showed that the 
concentration dependence of Tc in (Ga, Mn)As was well 
understood by the p-d exchange interaction if a correction 
to the local density approximation (LDA) is simulated by 
the LDA-I-C/ method with U = 4 eV4 

While these p-d exchange systems, in which the d- 
states of Mn impurities are practically localized, are well 
understood, there exist an even larger class of systems 
where the d-levels lie in the gap exhibiting impurity 
bands for sufficiently large concentrations. To these im- 
purity band systems belong (Ga, Mn)N, (Ga, Cr)N, (Ga, 
Cr)As, (Zn, Gr)Te, (Zn, Cr)Se and many others as shown 
by first-principles calculations^. Most of these systems 
are controversially discussed in the literature, and an un- 
ambiguous determination of the ferromagnetism has only 



been reported for (Zn, Cr)Te with a relatively high Gr 
concentration of 20% and a Gurie temperature of 300 
Kp2;. In particular in this class of materials, (Ga, Mn)N 
has been frequently mentioned as the most promising 
high-Tc DMS referring to the prediction of model cal- 
culations by Dietl et al,^ and ab-initio results by Sato et 
al.^ . Many groups have tried to fabricate ferromagnetic 
(Ga, Mn)N, but the experimental results are very con- 
troversial and confusing. After the first observation of 
the ferromagnetism of (Ga, Mn)?'^, many experiments 
followed, however reported Tc's are scattered between 
20 to 940 K,iojii2ii3ji4. Moreover, recently Ploog et al. 
observed spin-glass behavior in 7% Mn-doped GaN and 
suggested that the ferromagnetism observed in 14% Mn- 
doped GaN originated from Mn-rich clustersi^. Thus, 
the ferromagnetism in (Ga, Mn)N is still an open ques- 
tion which we reconsider in this letter. Ah initio calcu- 
lations by Akaii& and othero'^i^i^'^i^^i^Qi^" show that the 
magnetic properties of the above impurity band systems 
are dominated by double exchange mechanism and that 
the ferromagnetism is stabilized by the broadening of the 
impurity band. In the mean field approximation (MFA) 
high Tq values have been predicted [e.g., 350 K for (Ga, 
Mn)N with 5% of Mn, 500 K for (Ga, Cr)N with 5% of 
Gr, 400 K for (Zn, Gr)Te with 5% of Gr and so on) and 
the -^c-dependence of Tc on concentration c has been 
explained by band broadening^iii. Similar high, though 
slightly smaller, Tq values have also been obtained in the 
random phase approximation (RPA). 
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In this paper, we will show that a general obstacle 
for ferromagnetism exists in these dilute systems, in par- 
ticular in (Ga, Mn)N. Due to the large band gap the 
wave function of the impurity state in the gap is well lo- 
calized, leading to a strong, but short ranged exchange 
interaction, being dominated by the nearest neighbors. 
Therefore, for low concentrations, the percolation of a 
ferromagnetic cluster through the whole crystal cannot 
be achieved, so that a ferromagnetic alignment of the 
impurity moments cannot occur. Thus a paramagnetic 
or disordered, spin-glass hke, state is observed, in partic- 
ular for low concentrations. 

The electronic structure of DMS is calculated based 
on the local density approximation (LDA) by using the 
Korringa-Kohn-Rostoker (KKR) method. In this paper 
we focus on (Ga, Mn)N and (Ga, Mn)As as typical ex- 
amples for the double exchange and the p-d exchange 
systems, respectively. In these systems, Mn impurities 
distribute randomly at Ga sites in the host semiconduc- 
tor being described as (Gai_c, Mnc)X, where c is the Mn 
concentration and X refers to N or As. To describe the 
substitutional disorder, we use the coherent potential ap- 
proximation (CPA) . In this framework, all Mn impurities 
are equivalent and consequently we suppose a ferromag- 
netic alloy. It has already been shown that the mag- 
netic properties of metallic ferromagnetic alloys are well 
described within the CPA^^. While the CPA describes 
the electronic structure in the mean field approximation, 
we go beyond this approximation and explicitly calcu- 
late the exchange interaction Jy between two impurities 
at sites i and j, which are embedded in the ferromag- 
netic CPA medium. For the evaluation of we use the 
frozen potential approximations^ and apply a formula by 
Liechtenstein et al?^ . According to this formula, the to- 
tal energy change due to infinitesimal rotations of the 
two magnetic moments at site i and j is calculated us- 
ing the magnetic force theorem, and the total energy 
change is mapped on the (classical) Heisenberg model 
H = — Ei^j JijCiej, where ei is a unit vector parallel to 
the magnetic moment at site i, thus resulting in the ef- 
fective exchange coupling constant Jjj . This approach 
is already employed to estimate magnetic interactions 
in DMSs by Turek et al?^ and Bouzerar et alm^. For 
the present KKR-CPA calculations, we use the package 
MACHIKANEYAMA2000 coded by Akai2&. We assume 
mufhn-tin potentials and use the experimental lattice 
constants of the host semiconductors^. It has already 
been shown that the lattice relaxations in (Ga, Mn)N 
and (Ga, Mn)As are very small^"'^^'^^. Zinc blende struc- 
tures are assumed both for GaN and GaAs. In reality, 
GaN has Wurtzite structure. However, results for both 
structures are practically identical, because splitting of 
impurity bands due to symmetry lowering is small^*^ and 
disorder induced band width always overcomes the split- 
ting. The angular momenta are cut off at ^ = 2 in each 
mufhn-tin sphere. All calculations are performed for the 
neutral charge state of Mn, so that doping effects are not 
included. 
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FIG. 1: Calculated exchange interaction Jij in (a) (Ga, Mn)N 
and (b) (Ga, Mn)As as a function of distance. 



Fig. n shows the calculated exchange interactions Jy 
in (Ga, Mn)N and (Ga, Mn)As. As shown in the fig. [j- 
(a), in (Ga, Mn)N the interaction strength is strong, but 
the interaction range is short, so that the exchange cou- 
pling between nearest neighbors dominates. For exam- 
ple, nearest neighbor interaction Jqi in 1% Mn-doped 
GaN is about 13.5 mRy, while the other interactions are 
almost 2 orders of magnitude smaller than Jqi except 
for Jo4. Therefore, in this case the very large mean field 
value of Tq is mostly determined by Jqi. For higher 
concentrations, Jqi is suppressed and the interaction be- 
tween next nearest neighbors becomes negative, resulting 
a complicated structure in the distance dependence of 
the exchange interaction. Concerning to the mechanism 
of the ferromagnetism, it has already been pointed out 
that the double exchange mechanism dominates in (Ga, 
Mn)N where pronounced impurity bands appear in the 
gg^p4j7j9ji^^ It is intuitively understood that the exchange 
interaction in (Ga, Mn)N becomes short ranged due to 
the exponential decay of the impurity wave function in 
the gap. In contrast to (Ga, Mn)N, the exchange inter- 
action has long tails in (Ga, Mn)As in particular for low 
concentrations as shown in fig. n-(b). The qualitative 
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FIG. 2: Curie temperatures of nearest neighbor Heisenberg 
model in fee structure. Tq's are calculated by the mean field 
approximation (solid line), the random phase approximation 
(dotted line) and Monte Carlo simulation (crosses) . The per- 
colation threshold is 0.20 for the fee structure. 



difference in the interaction range between (Ga, Mn)N 
and (Ga, Mn)As is apparent from the figure. In (Ga, 
Mn)As the p-d exchange interaction becomes important 
as shown in ref4i. Since the extended hole state mediates 
the ferromagnetic interaction^, the interaction range is 
long ranged in p-d exchange systems essentially. Actu- 
ally the interaction extends farther than 3 lattice con- 
stants (20th shell). For higher concentrations, due to the 
screening of the pair interaction by the other impurities, 
interaction range becomes slightly shorter. 

As is well known the LDA predicts the position of lo- 
calized d-levels at too high energy. However, according 
to recent calculations by Shick et al.^^, the LDA-I-U cal- 
culations only slightly affect the impurity bands at the 
Fermi level in (Ga, Mn)N due to the extended nature of 
the anti-bonding t2 states of the impurity bands. There- 
fore, the LDA provides a fairly good description of the 
magnetic properties of (Ga, Mn)N. Even if the nearest 
neighbor interactions are changed in the LDA-I-U calcu- 
lations, this will not affect much the Curie temperatures 
for low concentrations, because only the longer ranged 
interactions are relevant due to the percolation effects. 
On the other hand, as we have already shown in ref.^, 
the LDA-f U calculations with U=4eV yield a different 
description of the magnetism in (Ga, Mn)As. This effect 
could change the calculated Tq values slightly, however, 
the exchange interaction in (Ga, Mn)As still remains long 
ranged and the basic argument of the following discussion 
is not affected. 

It is well known that the Curie temperature in 
the mean field approximation T^^^ is calculated as 
^B?c^^^ = ^cSi^oJoi, where fee is Boltzmann constant. 
As shown in this equation, evaluation of T^^^ does not 
require any information on the interaction range, because 
only the sum of the coupling constants appears in the 



equation. This simplification leads to significant errors 
in the calculated Tq of a dilute system with low concen- 
trations. This fact is easily understood by simple consid- 
eration and is known as the percolation problem^. Let 
us suppose a Heisenberg model with a ferromagnetic ex- 
change interaction only between nearest neighbors (near- 
est neighbor Heisenberg model), and consider what hap- 
pens when the system is diluted with non-magnetic sites. 
When the concentration of magnetic sites is 100%, we 
have a perfect ferromagnetic network. Due to the dilu- 
tion, the network is weakened, and for a concentration 
below a percolation threshold the ferromagnetism can- 
not spread all over the system leading to paramagnetic 
state since due to missing longer ranged interactions the 
moments can no longer align. Obviously this effect is not 
counted in the mean field equation for Tq, because the 
dilution effect is included only as a concentration fac- 
tor c in the equation. In case of the nearest neighbor 
Heisenberg model, the percolation thereshold Cp for fee 
structure is 20%^. In real cases the exchange interac- 
tion could reach beyond the nearest neighbors and the 
percolation threshold might be lower. However, in this 
paper we are interested in the concentration range well 
below the nearest neighbor threshold Cp. Therefore the 
exact Tc values could be much lower than the mean field 
values, in particular for the double exchange systems like 
(Ga, Mn)N where the exchange interaction is very short 
ranged (fig.UKa)). 

In order to take the percolation effect into account, we 
perform Monte Carlo simulations (MCS) for the effec- 
tive classical Heisenberg model. The thermal average of 
magnetization M and its powers are calculated by means 
of the Metropolis algorithm'^^. Due to the finite size of 
super cells used in the simulation, it is difficult to deter- 
mine Tc from the temperature dependence of {M{T)). 
In particular, when considering dilute systems, finite size 
effects and appropriate finite size scaling are of particu- 
lar importance for a correct and efficient evaluation of Tc 
by Monte Carlo simulations. To avoid this difficulty, we 
use the cumulant crossing method proposed by Binder—. 
This method uses the finite size scaling in the fourth order 
cumulant which is defined as C/4 = 1-{M^) /(3(M2)2). 
C/4 is calculated for various cell sizes and plotted as a 
function of temperature. If the cell size is larger than 
the correlation length, it can be shown that the Ui{T) 
curves for different sizes cross each other at three char- 
acteristic temperatures. Two of them are T = and 
T — 00, and the other is T = Tq. We use 3 cell sizes 
(6 X 6 X 6, 10 X 10 X 10 and 14 x 14 x 14 conventional fee 
cells) to carry out the cumulant crossing method for Tq 
calculations. For each temperature, we perform 240000 
Monte Carlo steps per site using every 20-th step for av- 
eraging. 

First, as a pedagogical example we show the calculated 
Tc for the dilute fee nearest neighbor Heisenberg model 
as calculated by MFA, RPA and MCS in fig.|21 For MCSs 
for dilute systems, we take 20 different random configu- 
rations of magnetic sites for the ensemble average. As 
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FIG. 3: Curie temperatures of (a) (Ga, Mn)N and (b) (Ga, 
Mn)As calculated by the MFA (solid lines), the RPA (dot- 
ted Unes) and the MCS (filled squares). For the MCS, the 
exchange interactions up to 15th shell are taken into account. 

shown in fig.|3 it is found that both MFA and RPA give 
reasonable estimations of Tc for c = 1 , with the RPA be- 
ing closer to exact MCS result. It has been analytically 
shown that for this model MFA gives upper limit of Tq 
and RPA gives lower limit'''^. However, for c < 0.7, MCS 
results are below RPA values and in particular below the 
percolation thereshold (cp = 0.20) the Curie tempera- 
ture vanishes: Tc — 0. Thus the serious deficiency of 
both MFA and RPA in the dilute concentration range is 
evident. 

Next, we show the calculated Tq values of (Ga, Mn)N 
(fig- 11(a)) and (Ga, Mn)As (fig. |l(b)) as obtained by 
the MCS from the Jy values in fig. ^ Thirty config- 
urations of Mn atoms are considered for averaging and 
Jij -interactions up to 15 shells are included; on the other 
hand, for the MFA and the RPA estimations interactions 
are included up to 72 shells. As shown in fig.|Sl-(a), very 
small Tc values are predicted for low concentrations in 



* Electronic address: |ksatoQcnip.sanken.os£ikarU.ac.jp| 



(Ga, Mn)N. MFA and RPA values are almost 2 orders of 
magnitude too large. Thus we find that the magnetism 
is strongly suppressed due to the missing percolation of 
the strong nearest neighbor interactions. Only the weak, 
longer ranged interactions satisfy the percolation require- 
ment, leading to small but finite Curie temperatures for 
5, 10, and 15 % of Mn. As shown in fig.lKb), due to the 
longer ranged interaction in (Ga, Mn)As, the reductions 
from the MFA are not very large, but still significant. 
Naturally these changes are larger for smaller concentra- 
tions. The Tc values of 103 K obtained for 5% Mn is in 
good agreement with the experimental values of 118 K 
reported by Edmonds et al2-. This value refers to mea- 
surements in thin films which are free of Mn-interstitials 
representing double donors. Including interactions be- 
yond the 15th shell, MCS could give slightly higher Tc 
values for low concentrations where the interactions do 
not converge within the 15th neighbors. At very high 
concentrations we expect our results to merge with the 
MFA and RPA values. 

In this Letter, we have shown by ab-initio calculations 
that (Ga, Mn)N shows no high-temperature ferromag- 
netism for low Mn concentrations. The strong ferromag- 
netic interaction of Mn nearest neighbor pairs does not 
become effective below the nearest neighbor percolation 
limit. The weak longer ranged interaction leads to a fer- 
romagnetic phase with very low Tc of several tens Kelvin. 
Therefore the experimentally observed very high Tc val- 
ues do not refer to a homogeneous ferromagnetic phase, 
but have to be attributed to small ferromagnetic MnN 
clusters and segregated MnN phases. Our results are 
of relevance for all DMS systems with impurity bands 
in the gap. To obtain higher Curie temperatures one 
needs longer ranged interactions and/or higher concen- 
trations. The latter requirement naturally points to II- VI 
semiconductors, having a large solubility for transition 
metal atoms. The observation of a Tc value of 300 K 
for (Zn, Cr)Te with 20% Cr^ is in line with these ar- 
guments. Similar results as presented above have been 
recently reported by a Swedish-Czech collaborationi24 
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